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Abstract. Multiscale transverse structures in the magnetic-field-aligned flows 
have been frequently observed in the auroral region by FAST and Freja satellites. 
A number of multiscale processes, such as broadband low-frequency oscillations 
and various cross-field transport effects are well correlated with these structures. 
To study these effects, we have used our three-dimensional multifluid model with 
multiscale transverse inhomogeneities in the initial velocity profile. Self-consistent 
generation of the low-frequency mode driven by local transverse gradients in the 
field-aligned ion flow and associated transport processes were simulated. Effects of 
particle interaction with the self-consistent time-dependent three-dimensional wave 
potential have been modeled using a distribution of test particles. For typical polar 
wind conditions it has been found that even large-scale (~ 50 — 100 km) transverse 
inhomogeneities in the flow can generate low-frequency oscillations that lead to 
significant flow modifications, cross-field particle diffusion, and other transport 
effects. It has also been shown that even small-amplitude (~ 10 — 20%) short-scale 
(~ 10 km) modulations of the original large-scale flow profile significantly increases 
low-frequency mode generation and associated cross-field transport, not only at the 
local spatial scales imposed by the modulations but also on global scales. Note 
that this wave-induced cross-field transport is not included in any of the global 
numerical models of the ionosphere, ionosphere-thermosphere, or ionosphere-polar 
wind. The simulation results indicate that the wave-induced cross-field transport 
not only affects the ion outflow rates but also leads to a significant broadening of 
particle phase-space distribution and tranverse particle diffusion. 


1. Introduction 

A number of in situ observations from OGO 5, Heos 

2, DE 2, and other spacecrafts indicate that large-scale 
transverse inhomogeneities in the field-aligned currents 
(FAC) and ion flows are ubiquitous in the ionosphere. 
The importance of the processes driven by transverse 
inhomogeneities in the ionospheric plasma flows has 
been shown both theoretically and from observations 
[D’Angelo et a/., 1974; Potemra et ah, 1978; Sugxura et 
at, 1984; Knudsen et ah, 1990; Loranc et al. , 1991; Ishii 
et ah, 1992; Horton et ah, 1993; Ganguli et ah, 1999], 

Recent high-resolution data from the Freja [Luhr et 
ah , 1994; Stasiewicz and Potemra , 1998; Knudsen and 
Wahlund, 1998] and FAST [McFadden et ah, 1998; Pe- 
na et ah, 1998] satellites showed that multiscale struc- 
tures in the parallel ion flows/beams and FAC with scale 
lengths from several hundred kilometers down to several 
hundred meters are a common feature of the auroral 


ionosphere. A number of effects, such as broadband 
low-frequency oscillations, transverse particle energiza- 
tion, and various multiscale cross-field transport pro- 
cesses, are well correlated with these highly structured 
parallel flows [Knudsen and Wahlund , 1998; Ergun et 
a/., 1998; McFadden et a/., 1998; Amatucci , 1999]. 

Adequate description of the multiscale processes as- 
sociated with inhomogeneous parallel flows usually re- 
quires three-dimensional (3-D) multimoment and mul- 
tifluid models. These models use much less computer 
resources than the more rigorous kinetic models and 
include important non-Maxwellian features which are 
not available in typical magnetohydrodynamic (MHD) 
global models. 

The one-dimensional (1-D) multimoment multifluid 
models have been successfully applied to a wide range 
of problems including large-scale modeling of the po- 
lar wind [Ganguli et al> 1987; Schunk and Sojka , 1989, 
1997], auroral return current regions [Ganguli and Pal- 
madesso , 1987], magnetosphere-ionosphere coupling 
[Mitchell et a/,, 1992], and potential structures genera- 
tion [S.B. Ganguli et of., 1994]. However, because of the 
1-D approximation, fluid instabilities due to transverse 
inhomogeneities in the parallel flows [e.g., D'Angelo , 
1965] and cross-field transport effects were ignored. 

The first 3-D fluid numerical simulations, which self- 
consistently addressed collective effects due to nonuni- 
form parallel flows and their influence on space plasma 
transport, were performed by Ganguli and Guzdar [1992], 
Ganguli [1995], and Ganguli et al [1999]. It was shown 
that the development of the shear driven instability 
[D'Angelo, 1965] can have a significant effect on iono- 
spheric plasma transport and saturated states. This ini- 
tial version of the model considered a uniform magnetic 
field and also assumed uniform and isotropic tempera- 
tures for all species. These studies were mainly focused 
on the dynamics of the hot magnetospheric plasma with 
application to the magnetosphere-ionosphere coupling 
problem. 

In this paper we have expanded on our previous work 
[Ganguli et al , 1999] to address low-frequency instabil- 
ities and associated transport effects due to transverse 
inhomogeneities in a field-aligned flow velocity for a typ- 
ical low-temperature polar wind plasma. Here we also 
demonstrate important wave-particle interaction effects 
through the use of a test particle distribution driven by 
a time-varying 3-D wave potential. The main focus of 
this paper is an analytical study and a 3-D simulation of 
the collective and transport effects due to the multiscale 
nature of the transverse structures in the parallel flow, 
as recently observed by the FAST and Freja satellites. 


For this purpose an initial parallel velocity profile with a 
typical large-scale transverse inhomogeneity modulated 
by small-scale perturbations of various amplitudes has 
been used. 

The paper is organized as follows. In section 2 we 
present a model description, the corresponding set of 
equations, and linear instability criteria. In section 3 we 
present our analytical and simulation results followed by 
discussions and conclusions in section 4. 


2. Model Description 

The model is designed to simulate the 3-D dynamics 
of a fully ionized (e~, H + ) plasma in a flux tube en- 
compassing a magnetic field line. The model is applied 
to the high-latitude topside ionosphere, and the simu- 
lation region extends in altitude from 1500 km to ~ 10 

Re- 

A detailed description of the model and its linear 
analysis was given by Ganguli et al [1999] and Mc- 
Carthy et al [1997] and is summarized as follows. The 
system of 3-D transport equations [Drake et a/., 1991] 
includes continuity and parallel momentum equations 
for the ions (equations (1) and (2)), an equation for 
the transverse velocity in the ExB drift approximation 
(equation (3)), and the relation between the electro 
static potential and plasma density (equation (4)) from 
the Boltzmann approximation for electrons along the 
magnetic field, 
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Here n and v are ion density and velocity, py = nuy, 
T e and T{ are electron and ion temperatures, m is ion 
mass, <f> is electrostatic potential ( E = — V<£), and B is 
the magnetic field. The initial equilibrium considered 
throughout this paper is described by a uniform density 
n = noi zero potential <t> = 0, and nonuniform paral- 
lel velocity u(x), where x is the transverse coordinate. 
This is an obvious stationary solution of (l)-(4). Equa- 
tions (l)-(4) are solved numerically using a trapezoidal 


time-stepping scheme, higher order spatial differencing 
and hyperviscosity [Guzdar et al , 1993; Ganguli et al , 
1999]. The boundary conditions are periodic in paral- 
lel z and transverse y directions. Constant values of 
all moments and zero gradients are maintained at the 
boundaries in the x direction. 

In this simplified version of the model the magnetic 
field is uniform and the temperature is uniform and 
isotropic. As discussed in section 4, these assumptions 
allow us to simplify the interpretation of the effects d ; s- 
cussed in this paper, while the essential results obtained 
here remain valid when more rigorous models axe used. 

In addition to the basic fluid model used earlier by 
Ganguli et al [1999], here we also solve 3-D equations 
of motion for a distribution of test particles (H + ions) 
moving in a 3-D time-dependent self-consistent wave 
potential. Full particle dynamics is considered along 
the magnetic field, while a gyrocenter approximation is 
used in the transverse direction. The gyrocenter ap- 
proximation saves a significant amount of CPU time 
since the ion gyroperiod is usually much smaller than a 
typical time step of fluid simulations. Since only long- 
wavelength low-frequency (much below the ion gyrofre- 
quency) modes are modeled here, the use of the gyro- 
center approximation in the transverse direction is an 
adequate approach to illustrate wave-particle interac- 
tions. The equations of motion used for the test parti- 
cles are given by 
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where X,Y, Z and V x ,V y , V z are Cartesian particle co- 
ordinates and velocity components, respectively. Com- 
ponents of the electrostatic potential gradient (V0) are 
obtained from the numerical solution of the fluid equa- 
tions (l)-(4) and linearly interpolated to a current par- 
ticle position X, Y, Z. Both sets of equations (l)-(4) 
and (5) and (6) are solved simultaneously. 

First, we present linear results for the D’Angelo insta- 
bility following D’Angelo [1965] and Catto et al [1973]. 
We consider equilibrium with a uniform magnetic field 
in the z direction, a uniform density, and a nonuni- 
form parallel flow with a gradient dv^/dx ^ 0 in the 
x direction. Considering perturbations proportioned to 
exp(— iut -I- ik v y + ik z z ) and linearizing (l)-(4), we ob- 
tain the following dispersion relation for the D’Angelo 
mode in the local approximation: 
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where cji = (a; — fc c D||o) is the mode frequency in the 
ion frame and c 8 = [T e + Ti/m ] 1 / 2 is the sound speed, 
and is the ion cyclotron frequency. A more rigorous 
version of (7) based on a kinetic formulation is given by 
G. Ganguli et al . [1994] and Gavrishchaka et al [1998, 
2000 ], 

For uniform parallel flow (dv^/dx = 0), (7) re- 
duces to the conventional dispersion relation for the ion- 
acoustic mode with the solution u\ = c s k z . However, 
when 



the solution of the dispersion relation (7) is an unstable 
mode with the growth rate 

7 = c,k z (A - 1) 1/2 , (9) 

and with a real frequency in the ion frame u)\ T — 0. The 
growth rate maximized over k z ( 7 m ) for a particular 
value of k y is readily obtained from (9) and is given by 
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Note that although we cannot adequately apply ex- 
pression (9) and/or use fluid simulations in a range 
{kyPi) > 1, we have already shown from a more rig- 
orous kinetic dispersion equation given by G. Ganguli 
et al [1994] and Gavrishchaka et al [1998] and by 2- 
D PIC simulations [Gavrishchaka et al , 2000] that the 
D’Angelo mode is unstable over a wide range of k y (from 
( kyPi ) < 1 to ( kypi ) » 1). However, as shown by 3-D 
fluid simulations [McCarthy et al , 1997; Ganguli et al , 
1999] and by our 2-D PIC simulations [Gavrishchaka et 
a/., 2000], the long- wavelength modes (k y pi <£ 1) are 
the main contributors to large scale transport effects. 
Therefore in the following we consider only effects as- 
sociated with long-wavelength modes. These are ade- 
quately described by our linear (7) and 3-D fluid simu- 
lation (l)-(4) models. 

In the next section, typical parameters of the polar 
wind plasma is used to estimate the instability growth 
rate from (10) for different velocity gradients consistent 
with observations. Instability development and asso- 
ciated transport effects have been studied using a 3-D 
numerical model based on equations (l)-(4) and (5) and 
( 6 ). 


3. Results 

We have investigated plasma collective effects and as- 
sociated transport phenomena due to transverse inho- 
mogeneities in the field-aligned ion flow. The following 


( 11 ) 


parallel velocity (u.) profile is used for this study 
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Here v z o is the magnitude of the velocity, Lt is the trans- 
verse scale length, A and A are a small-scale modulation 
amplitude and scale length, and the magnetic field is in 
the z direction. When A = 0, (11) describes a localized 
profile with a single transverse scale length ( L t ). This 
case will be used to model flows with relatively large- 
scale transverse inhomogeneities. When A > 0 and 
(A/2) <§: L t , (11) describes a filamented multiscale ve- 
locity profile, where the average transverse scale length 
is determined by Lt, while the local inhomogeneity scale 
length is determined by (A/2). The latter case will be 
used to model multiscale parallel flow profiles estab- 
lished by recent in situ observations [Luhr et a/., 1994; 
Stasiewicz and Potemra, 1998; Knudsen and Wahlund , 
1998; McFadden et a/., 1998; Ergun et a/., 1998; Pe- 
ria et a/., 1998]. Varying L t , A, and A, we can obtain 
profiles with different transverse scale lengths, scale sep- 
arations, degree of filamentation, and local velocity gra- 
dients dv z /dx. 

A profile example with L t = 60 km, v z0 = 20 km 
s” 1 , and A = 0 is shown by the solid line in Figure la. 
The same case but with A = 0.2 and (A/2) = 10 km is 
shown by the dashed line. The corresponding velocity 
gradient dv z /dx is plotted in Figure lb. It is clear that, 
although the average velocity gradient (determined by 
Lt) does not change when small-scale modulations are 
introduced, the local velocity gradients (determined by 
A and A) can be significantly increased even by a small 
amount of small-scale modulation. This is evident from 
(11), since dv z /dx includes two leading terras: one is 
oc 1/L e (average gradient) and the other is a A/ A (local 
gradient due to small-scale filamentation). Since the 
D’Angelo instability is local in nature and is driven by 
a local velocity gradient (dv z /dx), the modulation of the 
initial large-scale flow can have dramatic consequences 
as will be shown later. 

We begin our analysis by estimating the linear growth 
rates of the D’Angelo instability for a typical polar 
wind plasma with different transverse scale structures in 
the parallel velocity profile. The results are presented 
in Figure 2, where an (e-, H+) plasma is considered 
with T e = 1 eV, 77/+ = 0.2 eV, = 104 rad s“\ 
c s = 10.7 km s” 1 , and velocity magnitude v z o = 20 km 
s” 1 . In Figure 2a we plot the growth rate, optimized in 
k z , given by (10) for different values of the transverse 
localization Lt . We use the maximum value of dv z fdx 
as obtained from (11) for A = 0 (i.e., no modulation). 
Here k v = 1.14 x 10~ 4 m” 1 (i.e., A„ = 5.5 x 10 4 m), 


Figure 1 


Figure 2 




which will be used later in our 3-D simulations. One can 
see that even for large transverse scale sizes (~ 50-100 
km) the growth rate of the D’Angelo instability is sig- 
nificantly large. 

A more important feature is the growth rate sensi- 
tivity to the transverse inhomogeneity scale length. As 
seen from Figure 2a, a 2 orders of magnitude reduc- 
tion in the transverse scale length leads to an almost 2 
orders of magnitude increase in the. growth rate of the 
D’Angelo instability. Thus highly localized flow pro- 
files can be a very effective source of the low-frequency 
instability and associated transport effects. 

However, an isolated highly localized flow filament 
is not the most interesting flow configuration. First, 
it is not the most natural configuration encountered in 
the real space environment. Second, all collective and 
transport effects induced by this localized profile will 
occur in a small spatial region determined by the fil- 
ament transverse scale leading to minimal effects on 
the global plasma dynamics. A more natural config- 
uration confirmed by observations [e.g., Stasiewicz and 
Potemra , 1998; Peria et al , 1998] is a multiscale flow 
and/or FAC profile with multiple transverse structures 
of different scale. As shown later, this configuration is 
capable of inducing significant large-scale transport ef- 
fects in addition to localized processes. Equation (11) 
models such a profile with two different scales and mul- 
tiple filaments. In the more natural general case instead 
of one filament scale, one can use a statistical distribu- 
tion of different filament scales. However, to illustrate 
the main effect and for an easier theoretical interpreta- 
tion, (11) is an appropriate choice. 

Now we use (11) with v z0 = 20 km s -1 , L t ~ 60 km, 
and various modulation (filament) amplitudes (A), and 
scale lengths (A). The maximum local velocity gradients 
(dv z /dx) obtained from differentiation of (ll) are shown 
in Figure 2b as a function of A for (A/2) equal to 20, 10, 
5, and 1 km (solid, dashed, dotted-dashed, and dotted 
lines, respectively). We see a significant increase in the 
local velocity gradient when short-scale modulations, 
even with small amplitude, are introduced. 

Since the growth rate of the D’Angelo instability is 
determined by the local velocity gradient (shear fre- 
quency), the maximum growth rate of the instability 
(optimized over k z ) is readily obtained from (10) by us- 
ing dv z /dz from Figure 2b. The resulting growth rate 
for modulation scales used in Figure 2b is shown as a 
function of A in Figure 2c. As in Figure 2a, we use 
k u = 1.14 x 10“ 4 m" 1 . It is clear that even small- 
amplitude short-scale modulations can lead to a signif- 
icant increase in the instability growth rate. 


Short-scale modulation^ lead not only to the growth 
rate increase for a particular transverse wave number 
but also to a significant broadening of the unstable 
range in k y (or more generally k±) space. As was 
discussed by Gavrishchaka et ai [2000], the wide k y 
range of the simultaneously unstable modes leads to 
a broadband continuous wave spectrum when even a 
small transverse electric field ( E x ) is introduced. This 
effect is due to differences in frequency Doppler shifts 
(k v V e, where Ve is the ExB drift velocity) experienced 
by the modes with different k y . A transverse DC elec- 
tric field is often observed in space [e.g., Amatucci , 1999, 
and references therein]. The D’Angelo mode can be un- 
stable over a wide range of k y , especially in the presence 
of the short-scale filament at ion. Therefore the broad- 
band low-frequency wave activity observed in correla- 
tion with multiscale parallel flows and ion beams [e.g., 
Ergun et a/., 1998] can be at least in part explained by 
the described mechanism. 

We have shown that even small-amplitude short-scale 
modulation (filamentation) of the initial large-scale field- 
aligned flow can lead to a significant increase in the ef- 
fective growth rate of the D’Angelo instability. This is 
due to the increase in the local velocity gradient. How- 
ever, a linear calculation cannot show how small-scale 
filamentation of the initial steady state will change the 
overall instability development and associated trans- 
port effects. Linear analysis is also unable to clarify 
whether local enhancements of the instability and as- 
sociated transport processes will lead to a significant 
change in the global plasma behavior on scales large 
compared with the filament (modulation) scale. To an- 
swer all these questions, we use our 3-D fluid model 
(l)-(4) with a distribution of test particles moving in 
a time-dependent 3-D wave potential according to (5) 
and (6). 

We start our simulations of the polar wind (e-, H+) 
plasma with T{ = 0.2 eV, T t = 1 eV, and an initial flow 
profile given by (11). Initially, small density perturba- 
tions oc si n(k v y + k z z) are introduced to seed the insta- 
bility. The following flow parameters are used: v z o = 20 
km s’* 1 , Lt = 30 km, and (A/2) = 10 km. The mod- 
ulation amplitude (A) is varied to study the effects of 
short-scale filamentation. We use a set of wave num- 
bers k y = 1.14 x 10” 4 m" 1 and k z — 9.91 x 10" 8 m” 1 
(i.e., A v = 5.5 x 10 4 m and X z = 6.34 x 10 7 m). Even 
though single-mode seeding limits the number of un- 
stable modes, it allows us to clearly illustrate all the 
important effects and facilitates the unambiguous in- 
terpretation of the results. Simulations with random 
phase multiple mode seeding (not shown in this pa- 
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per) confirm all the main results presented here. In 
the following we compare instability development and 
associated transport effects induced by inhomogeneous 
parallel flow both with (A > 0) and without (A = 0) 
short-scale filamentation. 

In Figure 3a we present the time dependence of the 
wave potential amplitude for A = 0, 0.1, 0.2, and 0.3 
(solid, dashed, dot-dashed, and dotted lines, respec- 
tively). It is clear from this figure that even small- 
amplitude short-scale filamentation significantly accel- 
erates the plasma response to a structured field-aligned 
flow. Note, however, that the saturated state is rela- 
tively independent of A. 

As suggested by the linear analysis presented earlier, 
the main effects of short-scale modulation are due to 
the increase of the local velocity gradient leading to 
enhancement of the D’Angelo instability growth rate. 
This is confirmed by our simulations as shown in Fig- 
ures 3b and 3c. Figure 3b presents the maximum veloc- 
ity gradient of the initial profile (11) as a function of A. 
The growth rate obtained from the simulation curves 
shown in Figure 3a is given in Figure 3c by a solid line. 
This growth rate is equal to a 1/2 slope of the straight 
line that provides the best fit to the upper envelope of 
the curve ln^^^) in the linear regime. The dashed 
line represents the linear growth rate obtained from (9) 
using the dv z /dx values given in Figure 3b. It is clear 
that the A dependence of the growth rate obtained from 
simulations is close to that of the linear growth rate and 
is consistent with the A dependence of the maximum 
velocity gradient. The difference in the absolute values 
of the analytical linear growth rate and that obtained 
from the simulation may be due to artificial diffusion 
and viscosity terms used in the simulations for numer- 
ical stability. The artificial dissipation terms provide 
stabilization of short-scale fluctuations in 'the simula- 
tion system to compensate for the absence of natural 
(kinetic) dissipative effects in the fluid approximation 
[Guzdar et ai, 1993; Ganguli et a/., 1999; Gavrishchaka 
et a/., 1999]. The high-order numerical scheme used 
here [ Zalesak , 1981; Guzdar et al. } 1993] allows also to 
avoid excessive numerical diffusion typical for usual up- 
wind differencing methods [e.g., Press et a/., 1992]. 

In Figures 4a-4e the dynamics of the parallel veloc- 
ity profiles (averaged over y) for A = 0 (solid line) and 
A — 0.2 (dashed line) are compared. The profiles are 
given for t = 0 (Figure 4a), 1774 s (Figure 4b), 2129 s 
(Figure 4c), 2838 s (Figure 4d), and 3193 s (Figure 4e). 
We see that even small-amplitude initial short-scale fil- 
amentation of the profile leads to significant differences 
in the subsequent flow evolution due to the D’Angelo 
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instability. Figure 4f presents both profiles in the satu- 
rated state (solid line, t = 4258 s; dashed line, t = 3193 
s). This figure illustrates that saturated states for both 
initial loadings are quite similar. 

Figure 4 also clearly illustrates the stabilization mech- 
anism of the D’Angelo instability. Instability develop- 
ment broadens velocity profile decreasing velocity shear 
down to subcritical value which leads to saturation. In 
the absence of the explicit sources of dissipation, the 
spatially averaged velocity in the simulation box re- 
mains to be constant. On the other hand, there is a 
tangeable decrease in the average kinetic energy of the 
parallel flow due to wave generation. In the presented 
simulation the change in the average velocity between 
initial and saturated state was less than 1% (in both 
cases), while the change in the average kinetic energy 
was ~ 16% (without filamentation) and ~ 22% (with 
filamentation). 

So far we have presented only fluid diagnostics from 
our simulation. We find that in some respects the dy- 
namics of the test particles gives a better illustration 
of the effects induced by short-scale filamentation. In 
Figure 5 we compare the dynamics of the test particles’ 
phase space (V^, X) for the same parameters and time 
as in Figures 4a-4e. The cases with A = 0 and A = 0.2 
are shown in the left and right columns, respectively. 
Initially, all particles are located at the plane z = 0 
and have velocities corresponding to the local fluid ve- 
locity. All figures with particles presented in this pa- 
per illustrate the integrated picture of all particles from 
all z and/or y planes. One can clearly see that in the 
case of short-scale initial filamentation, particle diffu- 
sion in phase space begins significantly earlier and is 
much more intense compared to the case without short- 
scale modulation. We also see that at earlier times dif- 
fusion is localized and clearly correlates with the max- 
ima of the local velocity gradients (see Figure 5b) in 
accordance with the local nature of the D’Angelo insta- 
bility. However, at later times transverse mixing begins 
(Figure 5c) and, finally, the initial small-scale structure 
is completely destroyed and only large-scale transverse 
diffusion is evident (Figures 5d and 5e). Thus initial 
small-scale filamentation not only increases local cross- 
field transport but, more importantly, global transport. 

Similar effects are evident from test particle dynam- 
ics in configuration space (x,y) presented in Figure 6. 
As in Figure 5, cases with A = 0 and A = 0.2 are 
shown in left and right columns, respectively. The time 
for Figures 6b-6e is the same as for Figure 5b-5e. In 
Figure 6a we present a snapshot at t == 1419 s. For 
t = 0 there is a simple uniform distribution. Again, we 
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see that particle diffusion starts earlier and is signifi- 
cantly more effective in the case with initial short-scale 
filament ation. The local nature of the transverse diffu- 
sion is clear at early times (Figures 6a and 6b). Later 
diffusion becomes global and no sign of the initial small- 
scale structure remains (Figures 6c-6e). It is interesting 
to note that both phase space and configuration space 
evolutions of the test particle distribution (Figures 5 
and 6) are qualitatively very similar to those observed 
in a fully kinetic 2-D PIC simulation of the D’Angelo 
instability [Gavrishchaka et a/., 2000]. 

Finally, in Plates 1 and 2 we present the time evolu- 
tion of the 3-D velocity profile for A = 0 and A = 0.2, 
respectively. Color represents velocity in km s _1 . Snap- 
shots in Plate 1 are given at t =0 (Plate la), 3193 s 
(Plate lb), 3903 s (Plate 1c), 4967 s (Plate Id), and 
in Plate 2 at t = 0 (Plate 2a), 2129 s (Plate- 2b), 2838 
s (Plate 2c). In both cases we observe significant flow 
modifications as the instability develops. The maxi- 
mum flow velocity drops by ~ 35% (~ 50%) from its 
initial value for A = 0 (A = 0.2). Comparisons of 
Plates 1 and 2 (i.e., A = 0 and A = 0.2) reveal that not 
only does the initial filamentation accelerate instability 
development but it also leads to development of finer 
spatial structures, even though only single- wave per- 
turbation was used. The latter is mainly due to phase 
mixing of the modes excited in different local velocity 
gradient regions created by the initial short-scale mod- 
ulations in the velocity profile. 

4. Discussion and Conclusions 

We have presented results from the application of our 
3-D multimoment, multifluid model to a typical polar 
wind outflow with multiscale transverse structures in a 
field-aligned ion flow. The 3-D set of fluid equations al- 
lowed a self-consistent simulation of the low-frequency 
D’Angelo mode driven by the transverse gradient in the 
parallel velocity and its effects on transport processes 
in linear and nonlinear stages. We have also solved the 
equations of motion for a large ensemble of test parti- 
cles driven by the 3-D time-dependent wave potential 
obtained from the fluid model. This important addi- 
tion to a fluid model allowed us to demonstrate wave- 
particle interaction effects and particle transport due to 
low-frequency shear-driven waves. 

Our linear analysis and simulations show that for typ- 
ical polar wind conditions (T* = 0.2 eV, T e = 1 eV, 
Vjj = 20 km s" 1 ) even a large-scale (Lx ~ 100 km) 
transverse inhomogeneity in a parallel flow velocity can 
trigger the low-frequency D’Angelo instability with a 
typical growth time 1/7 ~ 3 min. This is a rather con- 
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servative estimate, since only the long-wavelength limit 
has been considered (A v > 50 km). The growth time 
is smaller for shorter wavelength. We also find that 
instability development leads to a significant modifica- 
tion in the flow (maximum velocity changes by > 30%) 
due to induced cross-field plasma transport. Analysis 
of test particle dynamics reveals that the growth of the 
D’Angelo mode also leads to substantial particle diffu- 
sion in phase space and configuratiqn space. The lat- 
ter results are qualitatively consistent with typical fea- 
tures observed in more rigorous but small-scale kinetic 
2-D PIC simulations of the D’Angelo instability [ Gavr - 
ishchaka et a/., 2000]. 

The main focus of this paper has been a study of the 
effects associated with multi-scale transverse structures 
in parallel flows. We found that even small- amplitude 
(10-20%) short-scale (~ 1-10 km) filamentation of the 
parallel flow leads to significant acceleration of the low- 
frequency D’Angelo instability (which depends on local 
velocity gradient values) and associated transport ef- 
fects. More importantly, we have shown that short-scale 
filamentation leads not only to local transport enhance- 
ment on the scales determined by filamentation but also 
significant transport intensification on the global scale. 

Our results are relevant to the interpretation of the 
recent high-resolution FAST and Freja satellite observa- 
tions of the transversely localized multiscale flows and 
currents correlated with broadband low-frequency wave 
activity and transverse transport effects [McFadden et 
al. } 1998; Ergun et a/., 1998; Knudsen and Wahlund , 
1998; Stasiewicz and Potemra , 1998; Peria et a/., 1998]. 
The reported results are especially topical since in many 
cases observed field-aligned currents are too small [e.g., 
Amatucci et a/., 1994; Amatucci , 1999, and references 
therein] to explain low-frequency wave generation by 
current-driven mechanisms [ Kindel and Ktnnel 7 1971]. 
On the other hand, we have shown here and recently 
by kinetic linear analysis and PIC simulations [ Gavr - 
ishchaka et al. } 2000] that the D’Angelo instability can 
lead to a broadband continuous wave spectrum and a 
rich set of associated transport effects without any field- 
aligned current. 

In the present study we used a simplified model with a 
uniform magnetic field and an isotropic plasma temper- 
ature to provide a simple and clear illustration of all the 
important effects. Preliminary results from the more 
complete model, which includes a dipole magnetic field 
and an anisotropic ion temperature, shows that for typi- 
cal parameters used in this paper the D’Angelo instabil- 
ity can be simultaneously operational over a wide range 
of altitudes (several Re) [Ganguli and Gavrishchaka , 



2000]. Therefore we expect that conclusions presented 
here will not change in any major way when the more 
rigorous model is applied. The complete set of the quan- 
titative results from the model with dipole magnetic 
field and anisotropic temperature [Gavrishchaka et a/., 
1999] will be presented in a future article. 

It should be also noted that in the presented sim- 
ulation quite conservative parameters have been used 
which resulted in the instability growth time 7" 1 up to 
10 min. As indicated by our linear analysis (e.g., Fig- 
ure 2), the growth time for typically observed plasma 
parameters can easily be in the range of several seconds 
or less. 

The instability analysis and simulation presented in 
this paper have been restricted to H-f plasma with one 
ion component. In reality, space plasma is multi-ion. 
For example, at lower altitudes it is (O-f, H-I-) plasma 
where 0+ is a dominant component. Although the de- 
tailed study of the D’Angelo instability in the multi- 
component plasma is beyond the scope of this paper, we 
can summarize our preliminary results on this subject 
obtained from a full kinetic dispersion relation [Gavr- 
ishchaka et a/., 1998, 2000] generalized to two-ion com- 
ponent plasma as follows. 

The effect of the addition of the O-f component to the 
H-f plasma critically depends on the shear in the par- 
allel velocity of the O-f component. For example, when 
O-f component is added to the H-f plasma and On- 
flow is shearless, then the growth rate of the D’Angelo 
instability due to the shear in the H+ flow slowly de- 
creases as 0+ concentration increases. However, even 
in a plasma with 80 - 90% of 0+ small shear in the 
H-f flow ( dv z fdx < can still excite the D’Angelo 
mode. On the other hand, when the shears in the O-f 
and H-f flows are comparable, even a small addition of 
O-f can lead to a dramatic increase in the instability 
growth rate, since the value {dv 2 /dx)/flo +y which en- 
ters dispersion relation, is much larger for O-f than that 
for H-f. Further research on this subject is warranted. 
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Figure 1. (a) Velocity profile u*(x) and (b) corre- 

sponding velocity gradient dv z /dx obtained from (11) 
for Lt = 60 km, (A/2) = 10 km, v z o — 20 km s” 1 , A = 0 
(solid line), and A = 0.2 (dashed line). 

Figure 1. (a) Velocity profile v z (x) and (b) corresponding velocity gradient dv z /dx obtained 

from (11) for L t = 60 km, (A/2) = 10 km, v z o = 20 km s” 1 , A = 0 (solid line), and A = 0.2 
(dashed line). 


Figure 2. (a) Linear growth rate, optimized in k z 

(equation (10)), for a typical polar- wind plasma as a 
function of transverse scale length in the parallel ve- 
locity profile. Here an (e-,H+) plasma is considered 
with T e = 1 eV, 77/4- = 0.2 eV, fl//+ = 104 rad 
s“ l , c a = 10.7 km s~\ k y ~ 1.14 x 10 -4 m -1 , and 
velocity magnitude v z q — 20 km s” 1 . The maximum 
value of dv z /dx obtained from (11) for A = 0 is used, 
(b) The maximum local velocity gradient ( dv z jdx ) ob- 
tained from differentiation of (11) as a function of A 
for L t = 60 km and (A/2) equal to 20, 10; 5, and 1 
km (solid, dashed, dotted-dashed, and dotted lines, re- 
spectively). Other parameters are same as in Figure 
2a. (c) The growth rate (optimized over k z ) obtained 
from (11) for modulation scales used (Figure 2b) as a 
function of A. The values of dv z /dx from Figure 2b are 
used. Other parameters are same as in Figure 2a and 
Figure 2b. 


Figure 2. (a) Linear growth rate, optimized in k z (equation (10)), for a typical polar wind 

plasma as a function of transverse scale length in the parallel velocity profile. Here an (e-,H+) 
plasma is considered with T e = 1 eV, T// + = 0.2 eV, f 1//+ = 104 rad s” 1 , c 8 = 10.7 km s” 1 , 
ky = 1.14 x 10“ 4 m" 1 , and velocity magnitude v z q = 20 km s" 1 . The maximum value of dv z /dx 
obtained from (11) for A = 0 is used, (b) The maximum local velocity gradient (dv z /dx) obtained 
from differentiation of (11) as a function of A for Lt = 60 km and (A/2) equal to 20, 10, 5, and 
1 km (solid, dashed, dotted-dashed, and dotted lines, respectively). Other parameters are same 
as in Figure 2a. (c) The growth rate (optimized over k z ) obtained from (11) for modulation 
scales used (Figure 2b) as a function of A. The values of dv z /dx from Figure 2b are used. Other 
parameters are same as in Figure 2a and Figure 2b. 


Figure 3. Results of the 3-D simulations of the polar 
wind (e-,H+) plasma with Tj = 0.2 eV, T e = 1 eV, 
and an initial flow profile given by (11). Initially, small 
density perturbations cc sin^y 4* k z z) are introduced 
to seed the instability. The following flow parameters 
are used: v z o — 20 km s _1 , L t = 30 km, and (A/2) = 10 
km. The modulation amplitude (A) is varied to study 
the effects of short-scale filamentation. One set of wave 
numbers k y = 1.14 x 10“ 4 m” 1 and k z = 9.91 x 10~ 8 
m -1 (i.e., A y = 5.5 x 10 4 m and \ z = 6.34 x 10 7 m) 
is used, (a) Time dependence of the maximum wave 
potential amplitude for A = 0, 0.1, 0.2, and 0.3 (solid, 
dashed, dot-dashed, and dotted lines, respectively), (b) 
Maximum velocity gradient of the initial profile (11) as 
a function of A. (c) Growth rate obtained from the 
simulation curves (solid line) shown in Figure 3a and 
the corresponding linear growth rate obtained from (9) 
for simulation parameters (dashed line) as a function of 
A. 


Figure 3. Results of the 3-D simulations of the polar wind (e-,H+) plasma with T* = 0.2 
eV, T e = 1 eV, and an initial flow profile given by (11). Initially, small density perturbations 
oc sin (k y y + k z z) are introduced to seed the instability. The following flow parameters axe used: 
v z0 = 20 km s’ 1 , L t = 30 km, and (A/2) = 10 km. The modulation amplitude (A) is varied to 
study the effects of short-scale filamentation. One set of wave numbers k v = 1.14 x 10’ 4 m” 1 
and k z = 9.91 x 10’ 8 m’ 1 (i.e., A y = 5.5 x 10 4 m and X z = 6.34 x 10 7 m) is used, (a) Time 
dependence of the maximum wave potential amplitude for A = 0, 0.1, 0.2, and 0.3 (solid, dashed, 
dot-dashed, and dotted lines, respectively), (b) Maximum velocity gradient of the initial profile 
(11) as a function of A. (c) Growth rate obtained from the simulation curves (solid line) shown in 
Figure 3a and the corresponding linear growth rate obtained from (9) for simulation parameters 
(dashed line) as a function of A. — 

Figure 4. Dynamics of the parallel velocity profiles 
(averaged over y) obtained from 3-D simulation for A = 

0 (solid line) and A = 0.2 (dashed line). The profiles 
are given for (a) t = 0, (b) 1774 s, (c) 2129 s, (d) 2838 s, 
and (e) 3193 s. (f) Both profiles are shown in saturated 
state (solid line - t = 4258 s, dashed line - t = 3193 s). 

Other simulation parameters are same as in Figure 3. 

Figure 4. Dynamics of the parallel velocity profiles (averaged over y) obtained from 3-D 
simulation for A = 0 (solid line) and A = 0.2 (dashed line). The profiles are given for (a) t = 0, 
(b) 1774 s, (c) 2129 s, (d) 2838 s, and (e) 3193 s. (f) Both profiles are shown in saturated state 
(solid line - t = 4258 s, dashed line - t = 3193 s). Other simulation parameters Eire same as in 
Figure 3. 

Figure 5. Comparison of the dynamics of the test 
particles’ phase space (V Z ,X) for the same simulation 
parameters and time as in Figures 4a-4e. The cases 
with A = 0 and A = 0.2 are shown in the left and right 
columns, respectively. 

Figure 5. Comparison of the dynamics of the test particles’ phase space (V Z1 X) for the same 
simulation parameters and time as in Figures 4a-4e. The cases with A = 0 and A = 0.2 are 
shown in the left and right columns, respectively. 

Figure 6. Test particle dynamics in configuration 
space (a:, y) for the same simulation as in Figures 4 and 
5. Cases with A = 0 and A = 0.2 are shown in left 
and right columns, respectively. The time for Figures 
6b-6e is the same as for Figures 5b-5e. In Figure 6a a 
snapshot at t = 1419 s is shown. For t ='0 there is a 
simple uniform distribution. 

Figure 6. Test particle dynamics in configuration space (x,t/) for the same simulation as in 
Figures 4 and 5. Cases with A = 0 and A = 0.2 are shown in left and right columns, respectively. 
The time for Figures 6b-6e is the same as for Figures 5b-5e. In Figure 6a a snapshot at t = 1419 
s is shown. For t — 0 there is a simple uniform distribution. 

Plate 1. Time evolution of the 3-D velocity profile 
obtained from simulation. Color represents velocity in 
km s’ 1 . Here A = 0 and other simulation parameters 
are same as in Figures 3-6. Snapshots are given at (a) 
t = 0, (b) 3193 s, (c) 3903 s, and (d) 4967 s. 

Plate 1. Time evolution of the 3-D velocity profile obtained from simulation. Color represents 
velocity in km s’ 1 . Here A = 0 and other simulation parameters are same as in Figures 3-6. 
Snapshots are given at (a) t = 0, (b) 3193 s, (c) 3903 s, and (d) 4967 s. 
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Plate 2. Time evolution of the 3-D velocity profile 
obtained from simulation. Color represents velocity in 
km s” 1 . Here A = 0.2 and other simulation parameters 
are same as in Plate 1. Snapshots are given at (a) t = 

0, (b) 2129 s, and (c) 2838 s. 

Plate 2. Time evolution of the 3-D velocity profile obtained from simulation. Color represents 
velocity in km s'" 1 . Here A = 0.2 and other simulation parameters are same as in Plate 1. 
Snapshots are given at (a) t = 0, (b) 2129 s, and (c) 2838 s. 
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Multiscale transverse structures in the magnetic-field- 
aligned flows have been frequently observed in the auroral 
region by FAST and Freja satellites. A number of mul- 
tiscale processes, such as broadband low-frequency oscilla- 
tions and various cross-field transport effects are well cor- 
related with these structures. To study these effects, we 
have used our three-dimensional multifluid model with mul- 
tiscale transverse inhomogeneities in the initial velocity pro- 
file. Self-consistent generation of the low-frequency mode 
driven by local transverse gradients in the field-aligned ion 
flow and associated transport processes were simulated. Ef- 
fects of particle interaction with the self-consistent time- 
dependent three-dimensional wave potential have been mod- 
eled using a distribution of test particles. For typical po- 
lar wind conditions it has been found that even large-scale 
(~ 50 — 100 km) transverse inhomogeneities in the flow can 
generate low-frequency oscillations that lead to significant 
flow modifications, cross-field particle diffusion, and other 
transport effects. It has also been shown that even small- 
amplitude 10 — 20%) short-scale (~ 10 km) modulations 
of the original large-scale flow profile significantly increases 
low-frequency mode generation and associated cross-field 
transport, not only at the local spatial scales imposed by 
the modulations but also on global scales. Note that this 
wave-induced cross-field transport is not included in any of 
the global numerical models of the ionosphere, ionosphere- 
thermosphere, or ionosphere-polar wind. The simulation re- 
sults indicate that the wave-induced cross-field transport not 
only affects the ion outflow rates but also leads to a sig- 
nificant broadening of particle phase-space distribution and 
tran verse particle diffusion. 
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